Abstract 



Basing on experimental data on irradiation-induced deformation 
of graphite we introduced a concept of diffuse domain structure de- 
vefoped in reactor graphite produced by extrusion. Such domains are 
considered as random continuous deviations of local graphite texture 
from the global one. We elucidate the origin of domain structure 
and estimate the size and the degree of orientational ordering of its 
domains. Using this concept we explain the well known radiation- 
induced size effect observed in reactor graphite. We also propose 
a method for converting the experimental data on shape-change of 
finite-size samples to bulk graphite. This method gives a more accu- 
rate evaluation of corresponding data used in estimations of reactor 
graphite components lifetime under irradiation. 
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1 Introduction 



Employment of graphite as neutron moderator as well as construction mate- 
rial has a long history starting with the first reactor created by Enrico Fermi. 
Artificial reactor grade graphites were designed for these purposes [Ij. The 
graphite components life time in active zone actually determines a life time 
of water-cooled carbon reactor as a whole, stimulating a study of radiation- 
induced evolution of reactor graphite properties. 

With respect to production technology, reactor grade graphites range 
from almost isotropic to substantially anisotropic. Here we imply both phys- 
ical and mechanical properties as well as irradiation effects. Widely used 
reactor grade graphites, such as Pile Grade A, ATR-2E, GR-280, produced 
by extrusion process exhibit a transversal-isotropic symmetry of physical and 
mechanical properties and radiation- induced effects [2H7j. Such symmetry 
requires five independent mechanical constants; in different representations 
these are elasticity constants, or compliances, or the set of Young's modulus, 
rigidity modulus, and Poisson coefficient [S]. 

Graphite rod subjected to irradiation undergoes transversely isotropic 
alternating-sign shape changes reaching several percent. The most important 
factors leading to non-uniform shape changes and hence, growth in internal 
stresses, are gradients of temperature and neutron exposure on the scale of a 
particular monolithic graphite element. Radiation- induced creep of graphite 
is not sufficient to compensate effects caused by shape changes and tempera- 
ture gradients. Finally, a complicated irradiation-induced evolution of elastic 
moduli along with shape changes may lead, under low creep, to a destruction 
of monolithic graphite components. 

In calculations of construction lifetime by using, e.g., the crack initiation 
and growth criterion, the determining factors are the data on evolution with 
neutron fluence of elastic moduli and relative shape changes in directions of 
an axis of symmetry (H) and in a plane of isotropy (-*-) [2H7l[9HT2]. The stan- 
dard way to obtain such data is by irradiation of finite-size graphite samples 
in materials testing reactors at a higher, compared with power reactors, neu- 
tron flux accompanied with periodical measurements of characteristics. In 
5-6 years the neutron fluence is acquired corresponding to 30-40 years lifetime 
graphite unit in operation conditions of power reactor. Such experiments are 
inevitably conducted under conditions different from operation ones, which 
has an inevitable influence on the data obtained. Necessary adjustments re- 
quire a development of the corresponding models which would help to convert 
the data obtained under higher neutron flux for finite-dimension samples to 
conditions of bulk graphite and power reactors. 
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2 Formulation of the problem 

In adjusting data on transversally isotropic reactor grade graphite irradiated 
in materials testing reactors one encounters three main problems, requiring 
a creation of models and a recalculation of corresponding data: 

I. Differences in neutron energy spectra. The problem is sufficiently 
simply solved by standard methods of neutron fluence recalculation 
to solid-state doze values in terms of displacements per atom (dpa) 
[5|fT3]. Note that such conversion is not universal. When heterogeneous 
nucleation of microstructural elements is important (for example, PWR 
pressure vessels) the conversion to dpa results in a partial loss of needed 
information. 

II. Noticeable differences in intensities of Prenkel pair genera- 
tion (dpa/ s). Morphology of reactor grade graphite has a compli- 
cated multiscale structure formed by the filler, which comprises the ob- 
jects with different degrees of order based on graphite micro-crystallites 
with the ideal crystal structure and dimensions of ~ 10^ nm, a binder 
with an isotropic micro crystalline structure, and an ensemble of micro- 
cracks [Il[SlinHI21[Ill[I5] ■ By the objects are meant formations of various 
scales from micro-crystallite complexes with various packing to larger 
formations including grains. As the dimension scale increases, the ob- 
jects become less anisotropic. An ensemble of microcracks is a result of 
relaxation processes at the scales with sufficiently high anisotropy level. 
At larger scales the anisotropy is too small to initiate the microcracks 
initiation. 

A driving force of radiation-induced effects in reactor grade graphite are 
processes taking place in an ideal crystalline structure of microcrystal- 
lite. Frenkel pair generation under the action of a neutron flux on crys- 
tal lattice results in the formation of a supersaturated two-component 
solid solution of interstitial atoms and vacancies. For a variety of rea- 
sons, the decay of quasi-2D supersaturated solid solution mainly leads 
to the formation of ensembles of interstitial dislocation loops of basal 
type, which, in turn, change microcrystallite shapes O [3 [121 US]. Be- 
low we list current understanding of the microstructural mechanisms 
of radiation-induced shape changes of microcrystallites: The formation 
of interstitial clusters, dislocation loops and new graphitic planes will 
cause an expansion of the microcrystallite in the c-axis [17\. Adjacent 
lattice vacancies collapse parallel to the layers and form sinks for other 
vacancies, causing shrinkage parallel to the graphite layers . Just 
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the change of crystalhte shape is a driving force for developing stresses, 
evolution of crack ensembles, and so on. At the bulk graphite scale, 
these effects are responsible for the shape changes, evolution of mechan- 
ical properties (partially) etc. Thus, juxtaposing radiation-induced ef- 
fects obtained under different generation rates of Frenkel pairs (dpa / s) 
requires solving the nontrivial problem on development of quasi-2D mi- 
crostructure in graphite crystal lattice with microcrystal self-consistent 
boundary conditions [T9] . 

III. Discrepancy between radiation-induced shape changes of bulk 
graphite and data obtained on finite-dimension samples. In our 

previous work [20j, based on numerous measurements of shape changes 
in cylindrical samples the two results were obtained: 

Initially circular cross-section of the sample under irradiation takes an 
increasingly pronounced elliptical shape and orientations of the ellipses 
vary randomly along the sample on a scale of 0.6 cm. The results of 
calculations of the relative volume changes using traditionally accepted 
methodology 

(Ay/y)y = (AL/L)" + 2(AL/L)^, (1) 

(where (AL/L)" and (AL/L)"*" are relative length changes of the sam- 
ples cut from the bulk material in directions parallel and perpendic- 
ular to the extrusion direction) increasingly diverge with irradiation 
doze from the results of the direct measurements of the relative vol- 
ume changes (see Fig. [1]). As it was shown [20] the reason for such a 
mismatch is that expression ([1]) only relates the relative elongations of 
samples to relative volume change while neglecting changes in cross- 
sectional areas. 

Basing on the presented experimental facts it was assumed that at macro- 
scale graphite can be considered as an transversally-isotropic elastic medium 
with irradiation-induced shape changes described by gradient tensor: 

/ 1 + (AL/L)^ ° \ 

Fm= [ 1 + (AL/L)^ . (2) 

V 1 + (AL/L)" / 

Here, (AL/L)"'"'" are the relative bulk graphite shape changes in parallel and 
perpendicular directions relative to the axis of symmetry (which coincides 
with the extrusion direction). But on a scale of ~ 0.6 cm, nonuniformities in 
shape changes are observed. 
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Figure 1: Relative volume change ^V/V (in %) as a function of fluence $ 
(neutrons / m^, En > 0.18 MeV) in samples irradiated at temperature 460 ± 
25° C. Two dash-dotted lines are experimentally measured volumes AL/L + 
2{Ad/d) in d = 6 mm diameter cylindrical samples of length L with both 
parallel (H) and perpendicular (^) cuts relative to the extrusion direction; 
single dotted line (T) shows the traditionally used combination ([1]), calculated 
using longitudinal experimental data for those two extrusions. For the sample 
with d = 8 mm, the result is qualitatively the same. 



Actually, to such regions without pronounced boundaries may correspond 
local (at a scale of / ~ 0.6 cm) random deviations of the axis of symmetry from 
the average direction in the bulk graphite. This representation was formalized 
by introducing a concept of "domain" - the region of size / possessing the 
gradient tensor of relative radiation-induced deformation in local coordinate 
system 

/ 1 + (Al/l)^ \ 

Fo = 1 + (A///)^ . (3) 

Here (A///)^ are irradiation-induced relative changes in sizes of free domain 
(i.e. mentially cut out of the elastic medium). The axis of symmetry of the 
domain is inclined at random angle 6 with respect to that of bulk graphite 
(see Fig. |2]). The domains are continuously transformed in space into each 
other on the scale of domain size I. It is impossible to uniquely identify the 
boundaries of such "diffuse" domains, which can be defined only statistically 
as regions in space with correlated directions of anisotropy. 

The longitudinal measurement data (AL/L)" , (AL/L)^ can be calculated 
by averaging over random angle 9 the projections of tensors Fq to the axis of 
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Figure 2: The random elastic medium with direction of anisotropy (sohd 
hne) continuously varying in space on the scale / can be visualized in several 
equivalent ways a) and b) as the sets of diffuse domains of size /. The axis 
e'g (anisotropy direction) of local coordinate system (e'^^, 62,63) related to 
domains is rotated by random angle 9 = 9 (x) with respect to the axis 63 
(orthotropy direction) of global coordinate system (61,62,63) related to the 
sample. 



symmetry of cylindrical samples. Such formalism leads us to writing a closed 
system of equations for obtaining (A///)|| , (Al/l)^ and sin^^, where 9 is the 
angular deviation of the domain axis of symmetry from the bulk graphite 
common axis of symmetry, by using experimentally measured parameters 
(AL/L)|, (AL/L)^, and e = ^Ef=i^fc where Ek = [c^max - (drain) k\/d is 
the flattening factor for k-th transversal cross-section of the sample. 

The developed model consistently explains the entire array of obtained 
experimental data. Taking into account that sample cross-sections are of 
the same scale [d = 0.6 cm and d = 0.8 cm) as domain size, the irradiation- 
induced change of domain shape in samples can be considered as almost 
free dilatation (parameters (A///)|| and (A///)j_). At the same time, if we 
consider the case of radiation-induced effects in bulk graphite, the presence of 
an ensemble of disoriented domains with changing shapes would give rise to 
internal stresses, which entails a conclusion that shape and volume changes 
in bulk graphite are not equivalent to those in finite-size cylindrical samples. 
Keeping in mind that we are aimed at determining shape changes of bulk 
graphite while experimental data refer to samples with cross-sections on the 
same scale as the domains size, we have to recalculate the experimental data. 
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This work is devoted to further development of model [20] in order to 
suggest a method for recalculating data obtained on small-size samples to 
radiation- induced shape changes of bulk graphite. In the model developed 
for transversal-isotropic medium we employ pure elastic approach rather than 
elastoplastic one as required for rigorous solution of the problem. This is 
explained by insufficient understanding of the mechanisms of graphite irra- 
diation creep related to the problem II above. Actually we obtain a lower 
bound for influence of disoriented domain structure on shape changes in bulk 
graphite. 

In section [3] we develop a new approach to the description of elasticity of 
such diffuse domain structure. We consider both the case of high anisotropy 
(see Fig. Ob) and uniform distribution of anisotropy directions (see Fig. Oc). 
Those cases correspond to different grades of reactor graphite [HHHB] ■ Macro- 




Figure 3: Domain structures in graphite: spherical shape of diffuse domains 
with transversely isotropic symmetry (shown by grid lines) and size / be- 
comes ellipsoidal (with sizes l± ^ due to irradiation-induced shape-change 
of free domains. Regular structure a) deforms affinely [L^ jL^ = l±/l\\). The 
domains in irregular domain structure b) and c) randomly deform due to 
interaction with their elastic environment, and the structure as the whole de- 
forms nonaffinelly {L-^/L^^ ^ i±/l\\)- In case b) of high anisotropy the macro- 
scopic deformation has the same transversely isotropic symmetry (L-*- ^ L"), 
but different amplitude compared with the domains. In case c) of isotropic 
distribution of domain directions the macroscopic deformation is isotropic 
(L-*- = L'I) even for anisotropic domain deformation ^ /y). 
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scopic description of diffuse domain structure on scales large with respect to 
domain size is developed in section HJ In section [S] we describe formation of 
diffuse domain structure in graphite produced by the extrusion process. We 
also estimate size and a degree of alignment of the emerging domains. 

3 Elasticity of diffuse domain structure 

In this section we introduce main concepts and formulate a physical model of 
elasticity of diffuse domain structure. We consider diffuse domains as more 
or less homogeneous regions with anisotropic elasticity. Such domains do not 
have sharp boundaries, and their direction of anisotropy continuously and 
randomly change in space on scale of domain size /, see Fig. [2l 

3.1 Strain of domain structure 

Deformation of solid can be described by strain tensor i, which is related to 
the gradient tensor as £ = {FF'^ — l)/2, where means transposition. Since 
the domain size is much larger than the scale of micro-cracks the strain tensor 
can be considered as a continuous function of coordinates x. We describe 
deformation of graphite by taking the state before irradiation as reference 
state. Thus defined strain tensor can be presented as the sum 5*°* = e'^ + e 
of two contributions: macroscopic deformation of sample 5^ (x) and local 
elastic deformation i (x) with components 

where components of elastic deformation vector. In equilibrium, this 

vector describes local deviations on the scale of domain size /, so that its 
average over the macro-scale vanishes. 

In the case of bulk graphite with free boundaries the macroscopic strain 
tensor is diagonal in global coordinate system related to sample (see 
Fig. [2]). In linear approximation in relative changes of sample sizes {AL/L)^ 
in principal directions a = 1, 2, 3 (see Fig. [3]b) it can be written as 

= - = i^L/L)^ 6.^ (5) 

Here and below, if possible, we shall drop the index M characterizing the 
bulk graphite. Gradient tensor Fm is defined in Eq. ([2]) and 6 is the unit 
Kronecker tensor. 

Graphite domains are randomly oriented in space. The local coordinate 
system related to the domains (see Fig. [2] b) can be obtained from global 
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coordinate system related to sample by means of rotations by three Euler 
angles Q: two polar angles if and ip and the azimuthal angle 6, characterizing 
the orientation of the local frame with respect to the graphite orthotropy 
axis. All three angles Q (x) = {ip (x) , 6 (x) , ip (x)) randomly vary in space 
with coordinates x on the scale of domain size /. Domain rotation by given 
angles fl = {ip, 9, ip) can be described by rotation matrix: 




(6) 



This matrix relates the strain tensors in global 
coordinate systems: 



e^^ + e and in local e 



loc 



-.tot 
'al3 



(7) 



Eq. ([7]) can be rewritten in matrix notations as 



.tot 






Here and below we use Greek indexes (a/3) to distinguish components of 
symmetric tensor 5^/3 = ^/3a and Latin indexes (i) when presenting this tensor 
in vector notation. The first three components of such vector determine 
normal strain Si = en (? 
vector determine shear strain, £4 = 
expressions for stress tensor are cTj 

0"13, 0-6 = '7l2- 

The rotation matrix R {Q) in Eq. 



1,2,3) and the other three components of this 
2ei3 and = 2ei2- Similar 
= 1,2,3) and 0-4 = 0-23,0-5 = 
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and where Vij are elements of the matrix (jH]). Below we use notations 
(Q) , (Q) and (Q) for inverse, transposed and inverse transposed 
rotation matrixes, respectively. 

3.2 Elastic energy of domain structure 

Energy of diffuse domain structure takes simplest form in local coordinate 
system, which continuously rotates in space (see Fig. [2]b): 

^ = IJ ^^E/^^- [^'"'^ (^) - . (^) - ^^]. (9) 

Irradiation-induced dilatation of freely dilating domains describes irre- 
versible deformation of the domain with free boundaries (which was cut out 
of the sample) due to irradiation. The dilatation is diagonal in local coordi- 
nate system related to the domain: 

e^^=(Fo) -5^p = {M/l)J^p (10) 

Here Fq is defined in Eq. ([3]) and (Al/l)^ are corresponding relative changes 
of free domain sizes in principal directions a = 1,2,3, see Fig. [3] a. The 
strain in local coordinate system i^"" (x) is related to the strain i in global 
coordinate system via Eqs. ([7]) and ([8]). Notice, that the above defined strain 
tensors have different origination: (Eq. [TOj) describes irreversible irradi- 
ation induced domain dilatation and creep, i is truly elastic deformation 
tensor and the macroscopic deformation e*^ (Eq. E]) has both inelastic and 
elastic contributions. 

Symmetric stiffness matrix c in Eq. ([9]) is nearly the same for all domains, 
but the axis of its symmetry randomly rotates with respect to direction of 
bulk graphite orthotropy axis. Components of local stiffness matrix c are 
elastic moduli of the domain. They change under irradiation and depend on 
neutron fluence and temperature (see section [H] below). Here we consider the 
case of transversely isotropic symmetry, when the matrix c in local coordinate 
system has the form 
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Substituting the local strain e''°^ from Eq. (|H]) to elastic energy (jH]) we 
finally get elastic energy of diffuse domain structure 

E = -l d^y Q, \r-' (e + i'') - i^] \r-' (i + - ^^1 (12) 

2 J L ' Jj L ' ij 

Using this elastic energy we describe below in section 14.41 three the most 
important cases of domain structures: 

Regular: Regular structure (see Fig. [3] a) deforms affinely with domain de- 
formation. In this case the equilibrium macroscopic deformation is the 
same as for free domains, Eq. ffTOjl . i^^ = e^"'^ = e^. 



Free: In the case of "unconnected" freely dilated domains = and 
the macroscopic deformation e*^ = Ri^ differs from deformation of 
individual domains only due to difference in domain orientations. 
Here the bar means averaging over domain orientation. 

Irregular: Domains of irregular structure (see Fig.|3]b) are in more cramped 
conditions with respect to the case of regular domain structure. Elastic 
environment tends to compress the domain with respect to the case of 
freely dilated domains and ^ . Equilibrium deformation i^'^ is 
determined by balance between elastic deformation of the domain and 
elastic response of its environment on this deformation. 

Uniform: In the case of uniform distribution of domain directions (see 
Fig. [3] c) the macroscopic deformation i^'^ is isotropic even for strongly 
anisotropic deformation of individual domains, e''"'^ ^ . 



4 Macroscopic description of diffuse domain 
structure 

In this section we develop macroscopic theory of reactor graphite elasticity. 
Although diffuse domains have different orientations, many such domains at 
spatial scales large with respect to domain size form an orthotropic and ho- 
mogeneous material from a marco mechanics standpoint. The main feature of 
diffuse domain structure is absence of cross-domain boundaries, near which 
the breakup of stress or strain may take place as in ordinary polycristalline 
solids. In the case of diffuse domain structure both stress and strain change 
continuously and the domains of diffuse domain structure directly feel the 
stress a (x) of the elastic enviroment. Under the influence of such stress the 
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randomly turned domains deform non-affinely with macroscopic deformation 
of bulk graphite (see Fig. [3] b and c) and the strain e of such irregular do- 
main structure gains large short wavelength (on the order of domain size /) 
components. 



4.1 Local stress 



Elastic energy f[T2l) can be expanded in powers of the strain e. We leave 
linear term in this expansion unchanged, and transform only its quadratic 
term using the equalitjfl 



~2 j d^^^.Sij(yi(rj + Jd^J2i 



(13) 



Here 

s{n) = R{n)siF (n) (i4) 

is local compliance matrix in global coordinate system and s = c'"'^ is local 
compliance matrix inverse to the local stiffness matrix c, Eq. (|TT|) . Let us 
show that if the stress a is found from extremum of the right hand side 
of Eq. (|T3|) . the latter will be equal to the left hand side of this equation. 
Substituting the solution of extremum equation for the local stress 

a, = ^ . Sr^' {Q) s, = Y,. \R-^ (^) cR~' (fi)! . . (15) 

3 3 -^13 

back to the right hand side of Eq. f[T3]) we really reproduce the expression 
for the left hand side of this equation. 

Using Eq. (|T3l) we can transform quadratic in e term in Eq. (|T2|) into 
corresponding linear term and rewrite elastic energy f|T2|) in terms of the 
local stress a (x) : 

E [a] = const — - d:si Sij [Q (x)] (x) aj (x) + 

2 J ^ij 

J d^J2, {a, (x) + af [n (x)] - a[ [fi (x)]) (x) (16) 



^Eq. is analogous to famous Hubbard-Strattonovich transformation, which is 

widely used to introduce weakly fluctuating order parameter (analog of the stress a (x)) 
in systems with strongly fluctuating dynamic variables (analoges of the strain e (x)). Ex- 
amples of such order parameters are magnetic moments in spin systems, superconducting 
order parameter in electronic systems and so on. 
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where is stress because of irradiation-induced deformation of the domain. 
In matrix notations 

af [n (x)] = i?-^ (^]) ce^ (17) 
Stress describes elastic response of medium on this deformation: 

a""' (fi) = R-^ {9) cR-' (n) i^' = (n) (is) 

Minimum of elastic energy f ITB]) gives standard equation of equilibrium in 
the presence of external stress (jj (x): 

Ea^K/^ + -5--a = o (19) 

To derive this equation we substituted elastic strain @ in Eq. (fT6|) . Variation 
of elastic energy for given stress tensor a takes the form 



5E = J d^J2^^ {a^, + - af,) 



d6uc 
dxR 



Integrating this expression by parts and equating volume term to zero, we 
find Eq. 

The solution of this linear equation can be written as a superposition of 
contributions from all domains. Each of these contributions slowly decays 
with the distance from the domain (as a power law at large distances). The 
domains have random orientations, and their total contribution to the stress 
is effectively self-averaged because of discussed above long-range character 
of elasticity. Therefore, in deriving equations of macroscopic elasticity we 
may consider the stress a (x) as weakly varying on the scale of domain size 
function. Since angles VL randomly alter in space, the constant stress a 
describes strain e (x) strongly varying in space on the scale of domain size, 
see Eq. ([T^. 

Note that Eq. f ITI?]) corresponds to minimum of energy fITB]) . Because of 
this minimum condition substitution of the approximate solution a (x) ~ 
const of stochastic equation (1191) in energy (|T6|) gives only minor corrections 
of second order in small-scale space variations of the function a (x). 



4.2 Derivation of macroscopic elasticity equations 

Below we will neglect the effect of weak spatial stress variations on the scale of 
domain size / on macroscopic elasticity of graphite. Than the integration over 
X in Eq. f|T6l) is equivalent to pre-averaging of the stress expansion coefficients 
Sij [Vt (x)] , 0"*^ [fi (x)] and a( [VL (x)] over characteristic spatial scale of the 
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function (t(x), and macroscopic free energy for slowly varying stress o" (x) 
takes the form 



EW] = -7; I d^y^, Sijai (x) aj (x) 

r 2J (20) 

+ j d^Y., (^) + (^) - (^)] (^) ' 

where bar means averaging over space. 

Equation of equilibrium (HM can also be transferred to its pre-averaged 
form with the aid of equality 



ax. w X — X — — — = — ax. — x 



dxfs J dx'p J dx'i^ 

^ ^dx'w(x'-x)F(x'^ ^^^"^ 



which is valid for any function F (x) and bell-like weight function w (x' — x). 
The weight function is centered at point x and has characteristic size larger 
than /. As the result of such averaging Eq. f|T9|) for slowly varying stress 
a (x) takes the form 

E ^ [^"/^ (^) + (^) - (x)] = (21) 

In equilibrium neighbouring domains can self-consistently adapt their de- 
formations since rigid directions of these domains are deployed by random 
angles to each other. Eq. [5T] should be solved with the boundary conditions 
on surface of bulk sample: 

E^ [^ap (x) + a'Jp (x)] = P« (22) 

where np are components of normal to the surface, and Pa are components 
of external force applied to the surface. This force can be considered as a 
substitute at the surface of corresponding volume source of deformation cr^^ 
in Eq. [HI 

Elastic energy fl20|) can be rewritten in standard form of macroscopic 
theory of elasticity by introducing "smoothed" on the domain scale / strain 
tensor e with components 

£.(x) = E.^^.^.(x) (23) 
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which is obtained by averaging the microscopic strain e with the above weight 
function w. 



Substituting Eq. (123|) into elastic energy (|20|) we can rewrite it in standard 
form for macroscopic elasticity: 



where C = S^^ is macroscopic stiffness matrix, af (x) and o"*^ (x) are av- 
erages of the corresponding values (1171) and ffTSjl . Minimum of this en- 
ergy reproduces equation of equilibrium of effective elastic medium fl?I]) . 
where the stress a is related to the strain e via Eq. f l2^ . Since we defined 
as equilibrium macroscopic strain (see Eq. [5]) in equilibrium we have 
aajs (x) = eai3 (x) = 0. Eq. (124|) with Ea/s (x) 7^ describing elastic response 
of graphite to external perturbation (elastic moduli, see next section). 

4.3 Elastic moduli of bulk sample 

The average compliance matrix S in Eq. fl20|) is found by space averaging the 
local compliance matrix 5* (fi), Eq. f|T4|) : 



Averaging of inverse stiffness matrix, and not of stiffness matrix itself, is 
related to the fact that domains deform non-affinely with the macroscopic 
strain, see Fig. [3]b [25]. Instead, they directly feel the applied macroscopic 
stress a (x) [26] . 

Space averaging in Eqs. fl20|) and fl25l) can be rewritten as averaging over 
all orientations of domains Q = {(p,6,ilj). Rotation angles ip and ip are 
randomly distributed in the interval (0, 2n) while probability distribution of 
the azimuthal angle 6 can be characterized by the second and fourth moments 
of sin 6: 



m2= sin^ 0^ (6) sin OdO, mi= sin^ {9) sin OdO (26) 





(24) 




(25) 
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The distribution function \I/ (9) of the angle 9 is calculated below in section [372] 
in case of graphite manufactured through the process of extrusion. In this 
case the fourth moment can be expressed in terms of the second moment 
1712, see Eq. HS] below. 

Calculating the angular average in Eq. (l25i) we get linear relation between 
elements of averaged and local compliance matrixes: 

Su = sii + I (-2sii + 2si3 + S44) 1712 + |s, 

5'l2 = Si2 + (-S12 + S13) 1712 + |s, 

5'i3 = si3 + i (sii + S12 - 3si3 + S33 - S44) i7i2 - is, (27) 

<S'33 = S33 + (2Si3 - 2S33 + S44) m2 + S, 

5'44 = S44 + I (6Sii - 2Si2 - 8Si3 + 4S33 - 5S44) 1712 " 2s, 

s = (sii - 2si3 + S33 - S44) m4 

In case of uniform distribution of domain directions (^^) = 1 substituting 
the moments 

ma = 2/3, = 8/15 (28) 

into the above equations we obtain moduli of macroscopically isotropic ma- 
terial (see Fig. [3]c): 

Sll = 5*33 = (8sii + 4si3 + 3S33 + 2S44) , 

5*12 = 5*13 = ^ (sii + 5si2 + 8si3 + S33 - S44) , (29) 
5*44 = 2 (Su - S12) = ^ (7sii - 5si2 - 4si3 + 2S33 + 5S44) 

Anisotropy and degree of disorientation increase with decreasing spatial 
scale and become large on scale of a stack of microcrystallites. Under the 
effect of thermal expansion of microcrystallites amplitude of stress on this 
scale can reach the level of crack initiation and growth. The characteristic 
scale of emerging microcracks is small compared to domain size. Therefore, 
relation (127|) between macroscopic and domain modules do not depend on 
evolution of microcrack subsystem and is determined only by the distribution 
of domain orientations. 



4.4 Free boundary conditions 

Solution of Eq. f l2T]) for samples subjected to external pressure P is a^'^ —a^ = 
(T^^* with diagonal tensor = P^ap- In case of free boundary conditions 
there is no stress at the boundary of the sample, a^^* = 0, and the equilibrium 
macroscopic strain e*^ is determined by the condition of zero average stress: 

= (30) 



17 



This condition corresponds to equilibrium between the average stress due to 
irradiation-induced deformation of domains, Eq. f lT7|) and the average stress 
due to response of elastic medium on this deformation, Eq. (fT7|) . Averaging 
the local stresses, Eqs. (|T7|) and (|T8l) . over Euler angles fl = {ip, 0, 9) we find 
linear matrix equation for the macroscopic strain (|5]): 



The solution of this equation is 

(AL/L)" = (A///)|| + A;" - (A///), 



(AL/L)^ = (A///)^ - (A///)^ - (A///)| 



(31) 



The factors k in this expression can be written in explicit form: 

/c" = [2 (4C44C|| — a) m2 + (c_L — C|[) 7777.2 

+2c|| (7-4044)7774]/;^, 

/C"^ = [2 (4C44C_L - A) 77l2 + 2 (c|j - C±) (7 + C_L - C||) 7^2 

+2cx (7-4c44)m4]/i^, 
i^: = 2A + [4c44 (c_L + 2c||) - 6 A] m2 - (c_l - cy)^ 7773 
+ (c± + 2c||) (7 - 4C44) 7714 

where 7772 and are moments of sin 9 (see Eq. and we use abbreviations 



(32) 
(33) 



C|| = Cii + C12 + Ci3, C_L = 2ci3 + C33, 

7 = cii + C12 - 2ci3, A = C33 (cii + C12) 



2r2 



Combining Eqs. (13T]) we find radiation-induced change of sample volume 
Ar/\/ = (AL/L)II + 2 (AL/L)^ 



= {M/l)^^ + 2 (A///)^ - kv {Al/l)^ - (Al/l) 
Amplitude of volume change is described by the factor ky- 

ky = k'^ — fc" = (c_L — C||) [8C447r72 — (37 + 2c_|_ — 2C||) 

+2 (7 -4044)7774]/ is: 



(34) 



rrir, 



(35) 



In case ky = {c± = cy) the relative volume change is the same as for 
medium with freely dilating domains 



{AV/V)p = (A///) II +2 (A///) 
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although such sample deforms strongly anisotropically (see Fig. Elb). In case 
kv 7^ the difference between those volume changes is related to deformation 
of domain shape due to elastic response of environment. The case ky > 
corresponds to compression of the domain by the medium, and ky < 
describes its dilation. 

Consider most important cases in more detail: 

Regular: Regular domain structure corresponds to case m2 = = 0, when 
the macroscopic strain coincides with the irradiation-induced strain, 

(AL/L)" = (A///)|| , (AL/L)^ = (Al/l)^ 

The bulk sample deforms affinely with domain deformation, see Fig. |3] 
a. 

Free: Stress in freely dilated domains is determined by the local condition 
(T^^ (x) = (x) that does not satisfy the equilibrium equation (HM due 
to long-range character of elasticity. The solution (x) = R[Q (x)] 
of the above condition takes into account only domain rotation and 
ignores elastic response of the medium on irradiation-induces deforma- 
tion of domains. Below for comparison we write down macroscopic 
deformation of such fictional medium of freely dilating objects: 

{AL/L)l = (A///)|| +m2 \{Al/l)^ - (A///), 

(AL/L)^ = (Al/l)^ - [(A//0^ - (A///),,] (36) 

and corresponding volume change of such medium: 



^ J F \ ^ J F \ ^ J F \^/|| \ ^ 

Irregular: At high anisotropy of domain directions (see Fig.[3]b) the factors 
k are proportional to the second moment 1112 ^ 1: 

fell ~ (4c44C||/A - 1) ma, k^ ~ (4c44C±/A - 1) ma (38) 

Nonaffine deformation of bulk sample is related both to random ro- 
tation of domains and elastic response of the medium on their shape- 
change. To separate both contributions compare Eqs. (13T]) and ( |38|) 
with Eqs. (15^ . The factors A;" and k^ are both proportional to ma, 
and similarly to Eq. fl36p they take into account the effect of domain 
disorientation on sample size. The difference between and ma is 

due to compression of domains by their elastic environment. 
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Uniform: In case of uniform distribution of domain directions (see Fig. E] 
c) with moments fl25]l we find from Eqs. fl52]) and fl55]) that the sample 
as the whole deforms isotropically under irradiation: 



{AL/iy = (AL/L) 
_fAl\ ^Jlc^ 




C_L + 2C| 



II V W II + 2c|, 



5 Formation of diffuse domain structure 

In this section we propose simple model of diffuse domains formation in 
graphite prepared by direct melt extrusion process. We also estimate the size 
of such domains and describe orientational alignment of domains in extrusion 
flow. 

5.1 Typical domain size / 

We will treat a mixture of pitch and coke particles at fabrication conditions 
as very viscous incompressible fluid. Velocity field v(x) of such fluid with 
density p can be described by Navier-Stokes equation: 



with boundary condition v = at channel walls. Here V is gradient, rj 
is dynamic viscosity of the fluid and force f initiates rotation of random 
graphite crystallites. Before extrusion the crystallites have random directions 
of anisotropy (with uniform distribution function (^) = 1); and they are 
rotated in the flow under influence of gradient of the force f . Mechanism 
of such rotation is different on spatial scales small and large with respect to 
characteristic correlation length 



where Re is effective Reynolds number: 

At small scales x < I main contribution to left hand side of Eq. (l39i) comes 
from viscous term —rjS/'^w ~ f . This term stabilizes the flow inside the 
correlation volume and describes regular laminar motion of the fluid 
on the scale /. Therefore, the rotation of anisotropy axes due to the 
force f is strongly correlated at small spatial scales x <l. 



p [d\/dt + (vV) v] - r/V V = f 



(39) 



/~ Re?7/ {pv) 



(40) 
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At large scales x > I main contribution to the turning force f comes from 
the term p(vV)v ~ f, which describes convective acceleration (insta- 
bility) of the flow and leads to randomization of motion at large dis- 
tances X > I due to amplification of the effect of random initial ori- 
entation of crystallites. This amplification also leads to a decrease in 
the Reynolds number Re with respect to case of ordinary structureless 
fluid. 

Using characteristic Reynolds number Re ^ 10 we find from Eq. fj40|) 
a simple estimate of domain size for viscous fluid with typical kinematic 
viscosity u = rj/p ~ 100 cSt (mm^ / s) and average velocity of flow v ~ 
10 cm / s: 

/ ~ Re iz/t; ~ 1 cm (41) 

This estimation is in good agreement with domain size / ~ 0.6 cm observed 
in GR-280 graphite [20j. Although our estimate of the Reynolds number Re 
needs additional experimental justification, Eq. f HT]) can be used to predict 
flow temperature and velocity gradient dependence of the domain size / in 
reactor graphites. 



5.2 Distribution of domain orientations 

Due to random shape and orientations of crystallites the turning force f in 
Eq. ( 139|) is also random, and coarse-grained at scale / domains in the flow 
experience rotations at random angles 66 during a small time interval 6t, 
see Fig. HI Such diffusion-like rotations can be described by orientational 
Smoluchowski equation for distribution function \E' (6) of domain angles 6: 



9^ (9) _ 1 d 

dt ~ ^^de 



D{e) sin0 



d^i (9) 



(42) 



with uniform initial distribution \E'o {6) = 1 before extrusion. In general, the 
angular diffusion coefficient D {6) = {66"^) /{26t) can depend on the angle 6. 
The drift term in this equation reflects the fact that the domain experiences 
an average torque r {6) due to the flow, which restricts its free rotation. Both 
opposite directions of domain anisotropy are equivalent, r {n — 9) = t (9) and 
therefore, the expansion of this function in a Fourier series has the form: 



' 'n>0 



Higher order corrections r„ with n > 1 are not important for our consider- 
ation and will be dropped below. Amplitude r = ri of the dimensionless 
torque is proportional to average strain rate of the flow. 
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Figure 4: Randomly turned coarse-grained at scale I domains rotate in the 
gradient flow v (x) at random angles 66 during a small time interval 6t and 
are oriented under the influence of orientation-dependent torque r {9). 



Steady state solution of Eq. fH21) at times t ^ D ^ has the form 

f7r/2 



cos^ ( 



sin ed9 (43) 



Substituting the distribution function fH5]) in Eq. fl2Bl) we flnd expressions 
for first two moments of this distribution: 

m2 = l-g' (r) , m^ = ml + g" (r) (44) 

In case of graphite obtained by extrusion at large velocity gradients the 
torque r is large, r 3> 1, and we get from Eq. (jUj) m2 ~ 1/r -C 1 and 
m4 ~ 2m2- In case of graphite formed under hydrostatic pressure average 
torque is small, r ^ 1, and we reproduce moments of uniform distribution, 
Eq. (l28l) . A typical value of the second moment of graphite GR-280 m2 — 0.3 
corresponds to intermediate case r ~ 4 between those two limits. 

In the Ref. [20] we show that the second moment 1712 can be directly ob- 
tained from analysis of experiments measuring irradiation-induced deforma- 
tion of small size graphite samples. The fourth moment can be expressed 
in terms of the second moment 1712 excluding parameter r from Eqs. fjUl) . 
Implicit dependence (jH]) of the fourth moment on the second moment can 
be interpolated for the whole range of torque r > by the following simple 
expression: 

2m'i , , 

rriA ~ -rrr (45) 

l + 3m2/2 ^ ' 
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6 Dose dependence 



Compliance moduli of a domain Sij can be estimated using experimentally 
measured macroscopic compliance moduli Sij. Difficulties in determining a 
set of five independent constants Sij characterizing elastic properties of trans- 
versely isotropic solid lie in the fact that, as a rule, only two Young's modules 
E\\, Ej_ and Poisson's ratio u are investigated in studies of irradiation effects 
of reactor graphite. In our case, we need to know all the five independent 
constants as well as their dose dependence. 

Experimental data for non-irradiated Pile Grade A reactor graphite [28] 
are shown in Table [H We expect that this graphite forms macroscopically 



Pile 
Grade A 


5*11 
2.150 


5'l2 

-0.127 


-0.123 


S33 
1.087 


3.333 


Mono- 
crystal 


0.011 


„m 

-0.00045 


„m 
■'IS 

-0.025 


„m 
■^33 

0.326 


„m 
^44 

4.35 


Domain 


2.5 


Sl2 

-0.12 


Sl3 

-0.19 


S33 
0.64 


S44 
2.73 



Table 1: Compliance moduli of Pile Grade A [28] and monocrystal non- 
irradiated graphites (in units of 10~^^ cm^/dyn). In last raw we show our 
estimate for domain complience moduli. 

homogeneous domain structure similarly to GR-280 graphite, since both of 
them are fabricated by the same extrusion method (see section O above). 
Values of moments m2 and 7714 and domain size I for these graphites may 
differ. 

Effect of irradiation on elasticity of Pile Grade A graphite has been stud- 
ied in details only for two Young's moduli E\\, E±_ and Poisson's ratios z/y and 
which are related to the macroscopic compliance moduli as [22|[25]: 

= 1/5*33, — 1/*S'44, l^W = — Si^/ S23, 

E^ = 1/Siu G^=l/[2 (Sn - ^12)] , = -3,2/ Sn. 

Dose dependence of entire set of compliance moduli has been established 
only for mono crystalline and pyrolytic graphites over a wide temperature 
range [21l[2SlEn]- It is noticed that only "small" elastic constants C33 and C44 
experience significant changes (increase in value) with irradiation dose. 

Detailed data processing for the elastic moduli £"11 and E± of ATR-2E 
graphite produced by extrusion, is given in Ref. [13]. It is shown that dose 
dependence of both Young's moduli i^y and E± is the same. The only differ- 
ence is in normalization factors: 

E||($) = i?,|/($), Ex ($) = Ex/ ($) (46) 
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Functional dependence of structure factor / ($) on neutron fluence $ is shown 
in Fig.|5]for different irradiation temperatures. Different experiments demon- 




10 20 30 40 10 20 30 




10 20 30 10 20 

<^ [dpa] <P [dpa] 

Figure 5: Dose $ (in units of displacements per atom) dependence / ($) = 
E ($) /E (0) of elastic moduli [13j of ATR-2E graphite at irradiation temper- 
atures 300°, 400°, 500° and 600°C (circles). Crosses show the approximation 
of the function / ($) by Eq. (HT]) . 



strate similar shape of dose dependence of the elastic moduli: a rather sharp 
increase, reaching maximum and subsequent decrease. Below we propose 
simple empirical expression for dose dependence of the structure factor / ($). 
We assume that there are two main contributions to the structure factor: 

/ ($) ^ 1 - ^9 (AV/V) + (47) 

due to evolution of subsystem of microcracks and due to irradiation hardening 
of crystallites. It is shown in Ref. [14j (Eq. 8) that the first of these contribu- 
tions to elastic moduli is proportional to relative volume change AV/V. This 
contribution is responsible for reversal of graphite shrinkage at high irradia- 
tion doses $ [HIEO]. Second contribution linearly grows with the irradiation 
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dose $ (different powers of $ give worse agreement with experiments). As 
one can see from Fig. O this simple dependence flTTI) describes adequately 
well all main features of dose dependence of the elastic moduli at different 
irradiation temperatures. Note also that the two fitting constants '(9 ~ 17 
and $* ~ 17 almost do not depend (at noise level) on irradiation tempera- 
ture. The universality of these constants is a further argument in support 
of equation fH7|) . Large value of dimensionless constant ^ is related to high 
contrast between elastic contributions of crystallites and microcracks. Simi- 
lar correlations between dimensional change and Young's modulus structure 
factor / (<l>) were also established for Gilsocarbon |27] . 

In order to determine dose dependence of all other elastic moduli, we first 
compare the Young's modulus G'||,± and z/|j _|_ for a typical Pile Grade 
A and single-crystal non-irradiated graphites, see Table [2J Due to random 



Pile 














Grade A 


9.2 


3.0 


4.6 


2.2 


0.06 


0.11 


Mono- 






eI^^ 








crystal 


909 


437 


30 


2.29 


0.04 


0.076 



Table 2: Young's {E) and shear (G) moduli of Pile Grade A ^28j and single- 
crystal non-irradiated graphites in units of GPa = 10^° dyn and correspond- 
ing Poisson's ratios (z/). Superscripts show preferential directions of crystal 
axes. 

orientation of crystallites main contribution to all macroscopic moduli Sij 
comes from "large" mono crystalline moduli and s^^ (see section 14.31 for 
description of orientational disorder effect). Macroscopic compliance moduli 
Sij are only weakly affected by "small" moduli and (see Table [1] 

for relative values of all moduli) since probability of finding basal planes 
along the direction of extrusion is small enough. In addition, anisotropy of 
elastic properties of reactor graphite is significantly smoothed out due to 
the presence of microcrack subsystem and binder [15]. The above described 
effects are responsible for observation of only moderate level of macroscopic 
orthotropy of graphite elasticity along with sharply pronounced anisotropy 
on crystallite scale. 

We conclude that dose dependence of the shear moduli of diffuse domain 
structure (see Fig.|2]) can be approximated by the same structure factor / ($) 
as for Young's modules fH6|) : 

G||(<l>) = Gx/(<l>), Gx($) = Gx/($) (48) 
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An important consequence of Eqs. ( 146 p and for the Poisson's ratios 

i^ll ($) = const, z/_|_ ($) = const 

is confirmed experimentally [2H]- Using Eqs. fH^ and f HS]) we find dose 
dependence of all macroscopic compliance moduli: 

S,, ($) = ($) 

where are corresponding values of those moduli in absence of irradiation, 
see Table [H 

Knowing experimental values of the moduli Sij Eqs. (127|) can be solved 
for the domain compliance moduli Sij 

s., ($) = s.,f-' ($) (49) 

Calculated values Sij in the absence of irradiation for Pile Grade A graphite 
are also presented in Table [TJ Substituting Eq. ( H^ into Eqs. (155]) 
and fl35p and using experimental value m2 — 0.3 and for GR-280 graphite j20| 
we estimate the factors 

A;ll ~ -0.17, ~ 0.18, ky ~ 0.35 (50) 

in Eqs. fl3Tl) and fl34l) for relative shape changes of bulk graphite samples. 
We conclude that bulk graphite sample deforms nonaffinally with domain 
deformation. 

Values of irradiation-induced dilatations (A///)|| and (A///)j_ can be es- 
timated measuring deformation of cylindrical samples (S) with diameter 8 
mm close to the domain size /. Relative changes in length of such samples 
with orientations coinciding (H) and orthogonal ("*") to the extrusion direction 
are [20]: 

(51) 

Here ml ~ 0.3 and ~ 0.7 are second moments of corresponding samples 
irradiated at temperature 550°C. 

Experimentally observed dose dependence of the relative length changes is 
plotted in Fig. El Using these data and substituting the solution (A//Z)|| j_ of 
Eqs. ( ET\i into Eq. (EI]) we also plot the predicted dose dependence of relative 
elongations of macroscopic samples. Note, that relative elongation of bulk 
samples is much higher than that of small size samples. In bulk samples the 
domains are in more crowded conditions in directions perpendicular to the 
extrusion direction (in these directions the deformation of domains due to 
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(AL/L)|^ = (A///)|| + ml'^ - (A///) 




Figure 6: Dose dependence of relative deformations (in %) of graphite irradi- 
ated at temperature T = 550° C in reactor BOR-60: The relative changes of 
the length of cylindric finite-size samples (S) and of macroscopic size samples 
(M) with orientations coinciding (||, solid lines) and orthogonal (_L, dashed 
lines) to the extrusion direction. 



orientational disorder is maximal). Therefore, at low doses, the relative strain 
in these directions is weaker, \AL/L\j^j < |AL/L|;^. Increase in logitudinal 
strain \AL/L\\.j > |AL/L|| is related to interaction between longitudinal 
and transverse elastic modes. These conditions can also be rewritten in a 
form {AL/L)j^ > (AL/L)^ and (AL/L)J^ < (AL/L)| that is also vahd for 
large doses. Since obtained results are sensitive to value of the moment m2, 
it is desirable to carry out further experiments to get better characterization 
of reactor grade graphites texture. 

7 Summary 

The main results of this paper are as follows: 

1. Basing on experimental data on radiation-induced deformation of 
graphite, we introduced the concept of diffuse domain structure, which is 
the reason of size effects in reactor graphite. In such structure there are no 
sharp boundaries between domains, and the direction of domain anisotropy 
undergoes random deviations from global orientation of the sample. We 
also derived equations of macroscopic elasticity, describing the deformation 
of solids due to irradiation-induced shape change of their domains having 
transversely isotropic structure. This theory can also be used to describe de- 
formation of stressed samples. A relation is obtained between elastic moduli 
of domains and macroscopic solid as the whole. We have proposed a scheme 
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for conversion of experimental data obtained on samples of finite size to de- 
scribe shape-change of bulk graphite. Such scheme is required for engineering 
calculations of graphite blocks integrity. 

2. We have presented simple model explaining origin of diffuse domains 
in graphite produced by the extrusion process. Our estimation of domain 
size is in good agreement with known experimental data [20j. This model 
can also be used to predict the dependence of the domain size on extrusion 
velocity and viscosity of visco-plastic medium. 

3. We also elaborated simple model describing orientational ordering of 
domains during extrusion. This model enables us to calculate the moments 
of domain's orientation, needed in description of elastic properties of reactor- 
grade graphite. 
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